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In long adaptive optics corrected exposures, exoplanet detections are currently 
limited by speckle noise originating from the telescope and instrument optics, and 
it is expected that such noise will also limit future high-contrast imaging instru- 
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Oh' ments for both ground and space-based telescopes. Previous theoretical analysis 

O ' have shown that the time intensity variations of a single speckle follows a modi- 

■ fied Rician. It is first demonstrated here that for a circular pupil this temporal 

intensity distribution also represents the speckle spatial intensity distribution at 
a fix separation from the point spread function center; this fact is demonstrated 
^ ■ using numerical simulations for coronagraphic and non-coronagraphic data. The 

^ ■ real statistical distribution of the noise needs to be taken into account explic- 

i/S ' itly when selecting a detection threshold appropriate for some desired confidence 

^ ■ level. In this paper, a technique is described to obtain the pixel intensity dis- 

^ ■ tribution of an image and its corresponding confidence level as a function of the 

^ . detection threshold. Using numerical simulations, it is shown that in the presence 

of speckles noise, a detection threshold up to three times higher is required to 
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obtain a confidence level equivalent to that at ha for Gaussian noise. The tech- 
nique is then tested using TRIDENT CFHT and angular differential imaging 
NIRI Gemini adaptive optics data. It is found that the angular differential imag- 
ing technique produces quasi- Gaussian residuals, a remarkable result compared 
to classical adaptive optic imaging. A power-law is finally derived to predict 
the 1 — 3 X 10~^ confidence level detection threshold when averaging a partially 
correlated non-Gaussian noise. 

Subject headings: stars: imaging, stars: low- mass, brown dwarfs, (stars:) plane- 
tary systems, instrumentation: high angular resolution, instrumentation: adap- 
tive optics, techniques: image processing, methods: data analysis, methods: sta- 
tistical 



Introduction 



Searchi ng for faint point sources around bright objects is a cha llenging endeavor. The 



atmosphere (IRoddiei 



1981 



ments optics ( Marois et all 



Racine et a] 



2003 



2005 



1999 



Macintosh et alll2005l ). telescope and instru- 



Hinkley et all 120071 ) produce speckles having a range 



of timescales that limit the direct detection of faint companions. From previous theoretical 
analysis, it is known that the speckle in t ensity temporal distribution is a modified Rician 
( lGoodmanlll968l : ISoummer &: Aimell2004l : iFitzgerald &: Graham! |2006| ). If a large number of 
uncorrelated speckle realizations are coadded, from the central limit theorem, then the final 
residual speckle noise follows a Gaussian intensity distribution. Since atmospheric turbulence 
produces random speckles that have a very short correlation time, a Gaussian distributed 
residual speckle noise is commonly assumed for ground-based adaptive optics (AO) long 
integrations and a detection threshold of 5cr is usually considered. 

However, careful residual noise analysis of AO images have demonstrated that long 
exposures are not limited b y random short-lived atmospheric speckles but by quas i -static 
speckles f lMarois et al I2OO3I : iMaroisI bood : iMasciadri eta] boosl : iMarois et a] boosl . boooh 
originating from the telescope and instruments. The speckle noise currently limiting high- 
contrast ground-based irnaging is thus very similar to that limiting space-based observations 
( ISchneider fc Silverstond 120031 ). The typical lif etime of ground-based quasi-static spe ckles 
has been found to be several minutes to hours (IMarois et all l2006l : iHinkley et all l2007l ) ; the 
noise in the combination of several images spanning ~1 hr is very similar to that in a single 
image (see Fig. [T]for an example; acquired with NIRI/Altair at the Gemini telescope). In 
this case, since the quasi-static speckle noise is well correlated for the entire sequence, the 
central limit theorem does not apply and the speckle noise in the final combined image will 
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be non-Gaussian. Sensitivity limits calculated assuming Gaussian statistics would have lower 
confidence levels (CL). Finding a robust technique to estimate proper sensitivity limits is 
fundamental to analyze adequately the sensitivity of an exoplanet survey as a function of 
angular separation. The contrast limit reached by a survey plays a central role in Monte 
Carlo simulati ons to derive exoplanet frequencies around stars and constrain planet forma - 
tion scenarios (jMetchevll2006l : ICarson et alll2006l : iKasper et all 120071 : iLafreniere et alll2007al ). 
Understanding the residual noise statistical distribution i s thus important fo r future ded- 
icated surveys of ne xt generation AO syst ems like NICI (IFtaclas et a 31200^, the Gemini 
Planet Imager fGPI. [Macintosh et all 120061 ) . the VLT SPHERE JPohlen et allbood ). and as 
well as future space observatories. 

In this paper, a new technique is presented to estimate sensitivity limits of a noise 
showing arbitrary statistics using a CL approach. The theory behind speckle statistics is 
summarized in § |2l Then § [3] presents a technique to derive detection thresholds using the 
probability density function and associated CLs. The technique is applied to simulated (§|1]) 
and observational (§ |5]) data to confirm the theory and to validate the technique. The effect 
of averaging a sequence of independent non-Gaussian noise realizations is discussed in § 16.11 
Concluding remarks follow in § [71 



2. Speckle Noise Statistics 



Following the work of lGoodmaru (119681 ) ; ISoummer fc Aimd (120041 ) ; [Fitzgerald fc Graham 



( I2OO6I ). the speckle intensity probability density function (PDF) for one location in the 
image plane and random temporal phase errors can be shown to be a modified Rician (MR) 
function. At a specific location in the image plane, the MR PDF Pmr{I) is a function of 
the local time-averaged static point spread function intensity Ic and random speckle noise 
intensities Ig'- 



where / is the point spread function (PSF) intensity (/ = Ic + h) and Io{x) is the zero-order 
modified Bessel function of the first kind. At a specific point of the PSF, if Ic ^ Is, relevant 
to Airy ring pinned speckles, the associated PDF is a Gaussian-like function showing a bright 
positive tail, while if Jc <C relevant to PSF dark rings or coronagaphic PSFs dominated 
by second order halo speckles, the noise distribution is exponential. The CL a for a given 
detection threshold d is simply obtained by: 



(2) 
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where p'^-^i is the mean-subtracted PDF. Fig. [2] illustrates the different possible regimes 
compared to a Gaussian intensity distribution. For a 5a detection threshold, where a is the 
standard deviation of the noise obtained using the robusLsigma IDL algorithmlll, a Gaussian 
distribution shows a 1 — 3 x 10~^ CL while a MR distribution show ~ 1 — 10"^ to 1 — 10"'^ 
CL. The MR distribution is thus producing much more false positive events. For example, 
consider a survey of many stars where each observation has a 500 x 500 A/D field of view 
(FOV, i.e. the 20" x 20" NIRI/Gemini FOV at H-band). If a 5a detection threshold is 
selected, the Gaussian noise distribution would lead to one false positive detection every 
four stars while the MR distribution would lead to ~250 to ~2,500 false positives per star. 
A detection threshold two to three times higher is required for the MR distribution to show 
the same CL as a 5cr Gaussian noise and the same number of false positive events. 

In the previous speckle PDF analysis, it was shown that the atmospheric speckle noise 
PDF is obtained by analyzing the temporal variation at one location of the PSF. For a quasi- 
static speckle noise, this approach is not adequate since the noise does not vary significantly 
with time. The quasi-static noise PDF can be derived using a very simple argument. If we 
consider a PSF produced by a circular aperture and if the PDF is obtained by analyzing 
pixels inside a narrow annulus centered on the PSF core, azimuthal quasi-static speckle noise 
variations Ig are produced with the same value of Ic (here, Ic is the unaberrated PSF and it 
is azimuthally symmetric for a circular aperture). The speckle noise inside a narrow annulus 
and from a single speckle noise realization thus shows the same PDF as a temporal speckle 
noise variation from random phase screens at any location inside the annulus. 



3. Experimental Derivation of the PDF and CL Curves 

A robust technique to derive sensitivity limits can be developed using CLs. The pixel 
PDF inside a specific region of the image is first obtained and the CL curve is then derived 
and extrapolated to estimate a local detection threshold. To avoid having too many false 
positive detections without missing possible faint companions, a CL of 1 — 3 x 10^^ {5a if 
Gaussian) is selected here. The basic steps to derive the PDF, to obtain the CL curves and 
to estimate the 1 — 3 x 10"'' CL detection threshold is summarized in tabled! 

The local PDF is obtained by producing an histogram of the pixel intensities inside a 
specific region of the image after subtraction of the mean pixel intensity over the region and 



^The robust_sigma algorithm uses the median absolute deviation as a first estimate of the standard 
deviation and then weight points using Tukey's Biweight; this algorithm provides another step of robustness 
to avoid biasing the standard deviation estimate if bad pixels are present. 
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division by the noise RMS of the image region. The CL curve as a function of detection 
threshold can be easily estimated by integrating the PDF inside the interval ±d (see Eq. |2]). 
Due to the limited number of resolution elements (one X/D for PSFs or a pixel for simulated 
noise images) in an image, the PDF and CL curves will be known up to a certain detection 
threshold. In theory, for the 1 — 3 x 10~^ CL detection threshold considered here, each area 
where the PDF needs to be estimated should have several million independent resolution 
elements. In practice, for images typically containing up to 500 x 500A/D (250,000 resolution 
elements), the PDF will be known only up to ~ 1 — 10^^ CL for Gaussian noise. A model 
fit using a analysis or a polynomial fit are required to extrapolate the CL curve and 
obtain the detection threshold corresponding to a 1 — 3 x 10^^ CL. Since the CL curves of 
various distributions are nearly linear in a semi-log(l — a) vs detection threshold plot (see 
Fig. [2]), we have chosen to use a polynomial fit due to its simplicity of implementation, its 
execution speed, and its accuracy. Due to non-hnear effects for detection thresholds near Oct, 
a linear fit is first performed for detection thresholds above 1.5a. If the detection threshold 
for a 1 — 3 X 10~^ CL is below 9a, a second order polynomial fit is used instead to better 
approximate the CL curve for quasi-Gaussian statistics. The CL extrapolation accuracy will 
be analyzed in the next section. 

4. Technique Validation with Simulated Data 

In this section, the PDF, the CL curve and the 1 — 3 x 10^^ detection threshold of 
simulated data are obtained. 



4.1. Simulated PDFs 

Simulated noise images using specific PDFs are used to test the algorithm in recovering 
the proper 1 — 3 x 10~^ detection threshold for known PDFs. To test the effect of the 
image area size on the CL extrapolation accuracy, images of various sizes are produced 
following a MR of Ic/h equal to 0.1, 1 and 10 (see Fig. [2j). For each size, 25 independent 
realizations are computed to derive the extrapolation accuracy. Fig. [3] and Table [2] show the 
CL extrapolation accuracy for simulated statistical distributions. In general, the algorithm 
slightly underestimates the 1 — 3 x 10^^ detection threshold for exponential statistics by 
~ 5%, but usually within the 2a error calculated for each area size. Typically, the bigger 
the area is, the more accurate is the detection threshold. To achieve a detection threshold 
accuracy of 10% for a 1 — 3 x 10~^ CL detection threshold, each PDF needs to be known 
up to a 1 — 10~^ CL (10,000 resolution elements per area). For instruments with smaller 
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FOVs (several hundred per several hundred X/D), an area of ~ 50 x 50X/D (2,500 resolution 
elements) would deliver a detection threshold accuracy of ~15% for all types of distribution. 
A solution to increase the detection accuracy of small FOVs would be to combine observations 
of several objects of similar magnitudes and observing conditions to increase the number of 
independant noise realizations in each area. 

4.2. Simulated PSFs 

The algorithm is now tested using simulated aberrated PSFs. For PSF observations, 
determination of the PDF is more complex. The speckle noise amplitude is decreasing 
with angular separation and the PDF may change with angular separation due to relative 
importance of random atmospheric speckles, photon, background and read noises. Since the 
speckle noise amplitude decreases with angular separation, a signal-to-noise ratio image is 
first obtained by dividing the pixel intensities, at each radius, by the standard deviation a of 
the noise at that radius (estimated using the IDL robusLsigma algorithm) . Finally, since CL 
are extrapolated, a compromise needs to be found between having a good radial sampling of 
the PDF and having sections of images big enough to adequately determine the PDF. 

PSF simulations are performed using Fast Fourier Transforms of complex 2048 x 2048 
pixel images with a 512 pixels diameter pupil; the full width half maximum (FWHM) of the 
PSF is 4 pixels. The pupil has uniform amplitude and includes A/160 RMS of phase errors 
generated using a power-law of index —2.6. The PSF images are then trimmed to 1024 x 1024 
pixels to avoid FFT aliasing effects. A non-aberrated reference PSF is subtracted to remove 
the Airy pattern and a signal-to-noise image is calculated. 

For simplicity, consider the calculation of the PDF within an annulus centered on the 
PSF. Since the presence of background or companion point sources inside that annulus could 
bias the statistics for real data (we are assuming that the background star density is such 
that only one or a few background objects are detected in the field of view (FOV) around 
any single target; cases with a high background star density will be discussed in section E^]), 
we have chosen to divide the annulus in three azimuthal sections containing 50,000 pixels 
each (~ 10,000 resolution elements, see Fig. |4]). The median PDF over the three azimuthal 
sections is calculated. Given the area of these sections, the PDF will be known down to a 
~ 1 — 10^^ CL for Gaussian statistics and the 1 — 3 x 10^^ detection threshold will be known 
to ~10% accuracy (see Tab. [2]). 

To further avoid cases where a point source is located at the border of two sections, the 
entire procedure is repeated by rotating the sections by 30 and 60 degrees which respect to 
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the PSF center, and the median PDF over the three orientations is finally obtained. This 
procedure is repeated at different angular separations. 

Since the PDF is estimated in large areas that may contain speckles with Jc ^ /s, /c ~ Is 
or Jc <^ Is, such technique returns an average PDF weighted by the various speckle noise 
contributions (pinned/unpinned speckles or Gaussian noises). Simulation with and without 
a coronagraph (simulated with a Gaussian pupil apodizer having a FWHM equal to a quarter 
of the pupil diameter) and for A/160, A/32 and A/16 RMS phase aberration are presented 
(see Fig. [5] and [6]). The algorithm clearly detects the MR distribution expected for a pinned 
speckle dominated PSF and the unpinned (exponential) speckle dominated coronagraphic 
PSF. 

The non-coronagraphic A/160 RMS simulations confirm that speckle noise follows a MR 
(required detection threshold of ~ lOo" for a 1— 3x 10~^ CL), as expected since pinned speckles 
are dominant for this case. As the quantity of aberrations increases, the ratio of pinned to 
non-pinned speckles decreases and the noise becomes exponential. For the Gaussian apodized 
case, since pinned speckles are strongly attenuated, the halo term dominates and the noise 
is more exponential. Note that none of these curves are expected to be fiat as a function 
of angular separation since the ratio of pinned to unpinned speckles is varying with angular 
separation, thus changing the pixel intensity distribution, and some noise is expected from 
the CL curve extrapolation (see Tab. [2]). Another simulation was performed using the A/160 
RMS case to show that if a constant Gaussian noise (background or read noise) is added to 
the image, the algorithm correctly detects the change of intensity distribution of the pixels 
at wide separations (see Fig. [7]). 

In high-contrast imaging observations, a partially correlated reference star PSF is usu- 
ally subtracted to remove a fraction of the quasi-static speckle noise. Such reference PSF can 
be obtained by observing a nearby target, b y acquiring the sam e star at another wavelength 



(simu ltaneous spectral differential imaging, iMarois et alll2005l ) or polarization (jPotter et al 



20011 ). or by building the reference usiri g images acquired with different field angles (angular 



differential imaging, iMarois et alll2006l ). Such a PSF subtraction is now simulated to esti- 
mate how its affect the PDF. The observed PSF / is simulated with a A/160 RMS phase 
aberration 0, with and without a Gaussian apodizer. The reference PSF Ij-^f is constructed 
by combining a perfectly correlated phase aberration a0, where a is a constant less than 1, 
with an uncorrelated part A0 such as: 

/=|FT(Ae'*)p (3) 
/,ef = |FT(Ae'("'^+^'^))|2, (4) 
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where the total noise RMS of cj) and acj) + A</) are equal and the ratio of the noise RMS of acf) 
and A0 is equal to 0.1, 1 and 10. Unless the background, read noise, random atmospheric 
speckles or photon noises are achieved by the reference PSF subtraction, the residual PDF 
is essentially unchanged for cases with and without a coronagraph (see Fig. [HD. 



5. Application to Observational Data 



The steps required to use the algorithm with observational data are similar to the ones 
described in Tab. H] and § 14. 2[ Only a few additional reduction steps are necessary. Fig. [9] 
illustrates the various steps of the technique using Gemini data. 

Beside the usual data reduction, deviant pixels, like diffraction from the secondary 
mirror support, must first be masked. Diffraction from the secondary mirror support usually 
produces a bright concentrated flux emanating from the PSF core along several azimuthal 
directions. Since this flux is not produced by quasi-static aberrations and is very localized 
in the image, if we include these pixels in the PDF, they will produced a bright positive tail 
in the PDF and the 1 — 3 x 10~^ CL detection threshold will be overestimated. 

Since the main science goal is to detect point sources, noise filtering is also applied to 
remove the noise that is not at the spatial scale of point sources. An 8x8 FWHM median 
filter is first subtracted from the image to reject large spatial period noises. Then, a 1 x 1 
FWHM median filter is applied to reject bad/hot pixels and smooth out the noise having 
spatial period below the resolution limit. The image is finally divided, at each radius, by 
the standard deviation a of the noise (again obtained with the IDL rohust_sigma algorithm) 
to obtain a signal-to-noise ratio image. 

The algorithm is first tes ted using data obtaine d at the Gemini telescope with the 
Altair adaptive opt ics system ( jSaddlemyer et all Il998h and the NIRI near- infrared camera 

fjHodapp et alll2000r ). These data are part of the Gemini Deep Pla net Survev f|Lafreriiere et all 

2007a ) that uses the angular differential imaging (ADI) technique f Marois et ai bood : iLafreniere et all 
2007bl ) to detect faint companions. This technique consists in acquiring a sequence of images 
with continuous FOV rotation. A reference PSF that does not contain any point sources 
is first obtained by combining images of the sequence, and the quasi-static speckle noise is 
then attenuated by subtracting the reference ADI PSF. The data for the star HD97334B 
(program GN-2005A-Q16), acquired on April 18, 2005 with good seeing conditions (Strehl of 
0.2 at H-band), are prese nted. These data h ave been reduced, registered and processed using 
the pipeline described in iMarois et all (120061 ) with the additional steps mentioned above, i.e. 
pixel masking and noise filtering & normalizing. Given that NIRI images are 1024 x 1024 



- 9 - 



pixels and PSFs have 3 pixels per FWHM, we have chosen the same areas as the simulated 
PSFs mentioned above (see § 14. 2p to calculate the PDF. For each region of the image, the 
pixel intensity histogram is obtained and then integrated to derive the CL curve. The CL 
is then extrapolated using a polynomial fit and the 1 — 3 x 10~^ CL detection threshold is 
estimated. The derived detection thresholds for a 1 — 3 x 10^^ CL are presented (see Fig. [TOl) 
for a single PSF, a PSF minus the ADI reference PSF, and the combined ADI-subtracted 
images. 

The algorithm is next tested using observational dat a obtained at the Canada-France- 
Hawaii telescope using the PUEO adaptiv e optics system (IRigaut et allll998l ) and TRIDENT 
near-infrared camera (IMarois et all l2005l ) . TRIDENT is a triple beam multi- wavelength 
(1.58, 1.625 and 1.68/xm with 1 % bandpass) imager built following t he simultaneous spectral 
differential imaging technique (IRacine et allll999l : iMarois et alll200Cll ). The technique consists 
in acquiring several images at different wavelengths and subtracting them to attenuate the 
speckle noise while retaining most of the flux of nearby companions. These data have been 
acquired as part of a direct imaging survey of stars confirmed to possess exoplanets from 
radial velocity analysis. The dataset of the star Ups And, acquired on November 14, 2002, 
is used here. Seeing conditions were relatively good and the Strehl ratio was of the order of 
0.5. The data have been reduced by subtracting a dark, dividing by a fiat field, correcting for 
bad/hot pixels, and registering the PSF at the image center. An optimized reference PSF was 
obtained using a star (Chi And) having a similar spectral type, magnitude and acquired at 
the same DEC and HA to minimize PSFs evolution from differential atmospheric refraction 
and flexure effects. Performances with and without the simultaneous reference PSF and the 
Chi And reference PSF subtractions are analyzed. Due to the limited FOV of TRIDENT, 
sections of 10,000 pixels (500 resolution elements, given that TRIDENT has 5 pixels per \/ D) 
are used to derive the PDF. Detection thresholds are thus known to ~15% (see Tab. [2]). 

It is clear that both TRIDENT and Gemini raw PSFs are limited by a non-Gaussian 
noise. For both TRIDENT/CFHT and NIRI/Gemini images, even after subtraction of a 
reference PSF, the residuals are still dominated by quasi-static speckles rather than averaged 
atmospheric speckle, background, read, or photon noises. Only the final combined ADI- 
subtracted image possessed a clear Gaussian-like noise. Fig.fTTlshows a visual example of a 5(T 
detection with and without a Gaussian distributed noise after introducing artificial 5cr point 
sources. For the Gaussian distributed noise, only the artificial point sources are detected with 
a detection threshold at 5cr0 while for the Gemini data (a MR distributed noise), numerous 



^Only approximately half of the artificial point sources are detected > since the artificial sources, 
being 5cr in intensity, vary in S/N by la RMS due to the underlying noise in the image. A ba detection 
threshold thus misses/detects ^ 50% of 5(t sources. 
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false positive sources are observed for the same detection threshold. If instead we select the 
1 — 3 X 10^^ CL detection threshold obtained by the technique described in this paper (here 
equal to lOa, see Fig. [10]), then only the artificial point sources are detected. It is interesting 
to note that the artificial source detection CL in the left and right panels of Fig. [TT] are the 
same. 



6. Discussions 

6.1. PDF Evolution with Quasi-Static Speckle Averaging 

It was shown in § [5] that the ADI technique produces a quasi-Gaussian noise. This is 
mainly due to the FOV rotation that occurs during the observing sequence; the residual 
noise is averaged incoherently when combining the images after FOV alignment. From the 
central limit theorem, it is thus expected that the noise in the final combined image shows 
quasi-Gaussian statistics. 

In this section, simulations are presented to estimate the number of independent speckle 
noise realizations required to converge to a quasi-Gaussian noise intensity distribution. Ran- 
dom noise images having 10^ resolution elements are created following an MR distribution 
having Id Is equal to 0.1, 1 and 10. The PDF and CL curves are calculated for a single 
realization up to the coaddition of 25 independent realizations (see Fig. [T2l) . Typically, ~20 
independent realizations are required for the MR distribution to converge, to ~20%, to a 
Gaussian distribution. Fig. [13] shows the detection threshold d for a 1 — 3 x 10^'' CL as a 
function of ngg, the number of independent noise realizations 

neff = n^^, (5) 

''"dcorr 

where n is the number of acquired images in the sequence, texp is the integration time per 
image and rdcorr the speckle noise decorrelation timescale (the equation is valid if Xdcorr > ^exp; 
if Tdcorr < 4xp then neff = n). These three curves can be well fit by a simple power-law of 
the form 

d{n) = [di - h]nj-'''' + 5, (6) 

where di is the detection threshold of a single image for a 1 — 3 x 10^^ CL. This equation 
is valid for all types of statistical distributions studied here. Eq. [S] can be used to predict 
the detection threshold required for a 1 — 3 x 10"^ CL and a statistical distribution with 
a known instantaneous PDF and speckle noise decorrelation timescale. If we consider the 
Gemini ADI observation {di ~ 13 for a single ADI-subtracted image, see Fig. [10]). for a 70 



minute observation sequence with rdcorr ~ 1-5 minutes (at 2" or 50A/Z^, see [Marois et al 
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20061 ). it is expected that the detection threshold for a 1 — 3 x 10 ^ CL of the final combined 



image [ries = 46.7) will be ~ 5.7 at 50X/D, in good agreement with the number derived 
with real images (~ 5.8, see Fig. [TUl) . 



6.2. Targets with a High Background Star Density 

In some cases, it is desirable to observe an interesting nearby target situated along the 
galactic plane or in front of the galactic bulge. In those areas, the high background stellar 
density implies that numerous background stars will be present in each of the areas defined 
to derive pixel intensity distributions. The detection threshold obtained will be affected by 
those stars since the algorithm would consider them as speckles. Several techniques can be 
used to remove the stars before estimating detection thresholds. A simple solution is to 
subtract these stars, using a non-saturated image of the primary, and mask any remaining 
contaminated areas with a not-a-number (NaN) mask. If the observations are obtained 
using the ADI technique then a star-free residual image can be obtained in the following 
manner. Prior to combining all the ADI-subtracted images together, instead of rotating 
them by the angle required to align their field of view, they are rotated by the negative 
of that angle such that all off-axis sources (the background stars) are eliminated by the 
median combination. As the amplitude of the rotation between the images is the same as 
for the "proper" combination, the effect of the median combination on the residual noise 
is expected to be the same and this star-free residual image can be used to estimate the 
noise distribution. Of course, the proper ADI combination of images must be used to search 
for companions. Another approach, still within the ADI framework, is to use the final ADI 
residual image to subtract the off-axis sources from each non-rotated ADI-subtracted image. 
Then these source-free images are rotated by the negative of the angle needed to align their 
field of view, such that their median combination eliminates the off-axis sources subtraction 
residuals. As for the previous technique, this source-free residual image should have the 
same residual noise distribution as the proper ADI residual image. 



7. Conclusion 

A robust technique was elaborated to estimate sensitivity limits using a CL approach. 
This technique correctly finds the expected MR intensity distributions of simulated and real 
PSFs, and properly detects a change of PDF as a function of angular separation. Experiments 
with simulated and observational data confirm the prediction of the theory that raw PSFs 
obtained with high-contrast imaging instruments are limited by a non-Gaussian noise. A 
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correction factor (up to 3) needs to be applied to detection limits found assuming Gaussian 
statistics to obtain the desired 1 — 3 x 10~^ CL detection threshold. Properly estimating 
this effect is important for future high-contrast imaging instruments for both ground- and 
space-based dedicated missions since a loss of a factor of three in contrast results in less 
sensitivity to low-mass exoplanets or, if a specific contrast needs to be achieved, integration 
times need to be at least nine times longer. It was shown that the ADI technique is the only 
observing strategy currently known that generates, intrinsically, a quasi-Gaussian noise at 
all separations where sufficient FOV rotation has occurred. A simulation has shown that it 
takes typically ~ 20 independent speckle noise realizations to produce an average speckle 
noise that shows quasi-Gaussian statistics. A general power-law is derived to predict the 
detection threshold required when averaging independent speckle noise realization of known 
PDFs and decorrelation timescale. 
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Table 1: Basic steps to derive the PDF, to estimate the CL curve and to determine the 
1 — 3 X 10~^ CL detection threshold. 



Step Action 

1 Define an image region 

2 Subtract the mean intensity 

3 Divide the pixel intensities by the noise RMS value 

4 Obtain the pixel intensity histogram (PDF) 

5 Integrate PDF between ±d to obtain the corresponding CL curve 

6 Perform a polynomial fit on the CL curve 

7 Extrapolate the CL curve to derive the 1 — 3 x 10~^ CL detection threshold 
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Table 2: The 1 — 3 x 10^^ CL detection threshold d derived by extrapolating different type 
of MR statistics as a function of the number of resolution elements. MROl is a modified 
Rician distribution with Ic/h — 0.1 while MRl and MRIO are for Ic/h — 1 and h/ L — 10 
respectively. Both the mean and standard deviation of the detection threshold are obtained 
by analyzing 25 independent noise realizations. 



Statistics 


Number of 


Expected d 


< ci > 


Stddev((i) 




res. element 


(^) 


(^) 


('^) 


Gaussian 


10^ 


5.0 


5.4 


1.1 




10^ 




5.33 


0.40 




10^ 




5.06 


0.29 




10^ 




5.06 


0.11 


MRIO 


10^ 


7.7 


9.3 


1.9 




10^ 




8.20 


0.95 




10^ 




8.02 


0.59 




10^ 




8.02 


0.66 


MRl 


10^ 


13.5 


14.9 


2.5 




10^ 




13.95 


0.98 




10-5 




13.44 


0.50 




10^ 




13.48 


0.37 


MROl 


10^ 


18.2 


18.7 


2.6 




10^ 




17.2 


1.4 




10^ 




17.13 


0.40 




10^ 




17.18 


0.41 



Fig. 1. — Gemini NIRI/Altair quasi-static PSF of the star HD97334B. A single 30s exposure 
is shown (left) as well as the median of 90 30s exposures (right). The inner saturated region 
has been masked. A symmetric radial profile has been subtracted to highlight the speckle 
noise. FOV is 20" x 20". Both images have the same intensity range and linear gray scale. 
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Fig. 2. — Left panel: PDF for a Gaussian distribution (solid line) and modified Rician with 
Ic/Is — 10 (dashed line), Ic/h — 1 (dot-dashed fine), and Ic/h = 0.1 (dotted fine). Right: 
corresponding CL as a function of detection threshold. 
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Fig. 3. — The PDF (left panel) and CL curves (right panel) obtained from a simulated noise 
image (10® resolution elements) generated with an MR intensity distribution compared to 
the analytical PDF and CL curves. In the left panel, the derived PDFs for Ic/h = 0.1 (thick 
line), Ic/Is — 1 (thick dot-dashed line), and Ic/h — 10 (thick dashed hne) are shown. The 
three thin hnes are the expected PDFs for the three same cases. The curves for the right 
panel are the same. The three thin dotted lines in the right panel are the extrapolated CL 
curves. 




Fig. 4. — Selected grid to calculate PDFs. 
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Fig. 5. — PDFs and CLs for simulated PSFs. The dotted line is for a Gaussian noise dis- 
tribution while the dashed and solid lines are respectively for simulated aberrated PSFs 
without/with a Gaussian pupil apodizer. The right panel shows the corresponding GL in- 
cluding the power-law fit (thin dotted lines) to estimate the 1 — 3 x 10~^ GL detection 
threshold. 
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Fig. 6. — Required detection threshold to obtain a 1 — 3 x 10^^ CL for simulated aberrated 
PSFs without (left) and with (right) a Gaussian pupil apodizer. The solid, dashed and dotted 
lines are for A/16, A/32 and A/160 RMS of phase aberration (generated with a power-law of 
index —2.6), respectively. 
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Fig. 7. — See Fig. [6l A random Gaussian pixel-to-pixel noise is added to the PSF image 
such that it dominates at wide separations. 
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Fig. 8. — See Fig. [61 A partially correlated reference PSF is subtracted (see text for details). 
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Fig. 9. — The different steps involved in estimating a detection tfiresliold in a specific region 
of tfie PSF using tfie confidence level approach. Panel 1 shows a Gemini CH4-short saturated 
PSF image that has been reduced and registered to the image center. Panels 2, 3 and 4 
show respectively the same PSF image but with the secondary support structure diffraction 
masked, after noise filtering, and after noise normalizing. Panels 5 and 6 show respectively 
the pixel intensity distribution (PDF) and corresponding CL curve inside a typical region of 
the PSF shown in panel 4. An extrapolation of the CL curve gives the 1 — 3 x 10~^ confidence 
level detection threshold for that region (here approximately 12.6a). Steps illustrated by 
panels 5 and 6 are repeated for all regions of the PSF and after rotating the regions by 30 
and 60 degrees to ehminate the bias resulting from point sources located at the edge of two 
regions. 
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Fig. 10.— Left panel: TRIDENT SSDI required detection threshold for a 1 - 3 x 10"^ 
CL of the star Ups And acquired at the Canada-Prance-Hawaii telescope. Solid, dashed 
and dotted lines show the required detection limit for respectively the TRIDENT PSF, the 
SSDI simple difference subtraction (subtraction of two images acquired simultaneously at 
two wavelengths) and the SSDI simple difference minus the reference SSDI simple difference 
image of a reference star. Right panel: same plot but for the ADI technique and the star 
IID97334B acquired at Gemini. Solid, dashed and dotted lines are for respectively a sin- 
gle PSF image minus a symmetric radial profile (spiders masked), a single ADI reference 
subtracted image and the total ADI median combined image (median of 90 ADI-subtracted 
images) . 
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Fig. 11. — 1 — 3 X 10"'^ confidence level detections for a Gaussian distributed noise and 
Gemini data. After analyzing the noise statistical distribution, artificial point sources are 
added at 0.9, 1.6 and 2.5" from the image center at PA of 20, 110, 200 and 290 degrees from 
the vertical axis and each having an intensity equal to 5a of the noise. Panel lA shows the 
image generated using a Gaussian noise (display with a linear intensity range between -la 
and 5a). Panel IB shows the pixels that are higher than 5a. Panels 2A and 2B are similar 
but for the Gemini PSF. Panels 3A and 3B are again similar but for artificial point sources 
now having an intensity equal to lOcr of the noise, the detection threshold derived using the 
CL technique. FOV is 5.65" x 5.65". Panel 3A is displayed with a linear intensity range 
between -2a and lOcr. 
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Fig. 12. — Extrapolated CL curves (&om 10 resolution element simulations) as a function of 
the detection threshold for various Ues- Upper left panel is for an MR statistical distribution 
with Ic/Is = 0.1 while the upper right and bottom panels are for Ic/Is = 1 and Ic/Is = 10 
respectively. The solid thick line is for a single noise realization while the thick dashed, 
dot-dashed and dotted lines are for the average of 2, 5 and 25 independent noise realizations. 
The thin dotted line is the CL curve for a Gaussian noise intensity distribution. 



-29- 




Fig. 13. — The 1 — 3 x 10~" CL detection threshold as a function of Ues, the number of 
independent noise reahzations. The thick hne is for an MR statistical distribution with 
h/h = 0.1 while the thick dashed and dot-dashed hues are for Ic/Ig — 1 and Ic/Is — 10 
respectively. The thin sohd hne is for Gaussian statistics. The three thin dotted hues are 
the predicted detection thresholds derived from the fitted power-law (see text for details) . 



